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We search for the decays B°_-^ p°p°, B° -> p°/o(980), and 5° -> / (980)/ (980) in a sample of 
about 384 million T(4S') — > BB decays collected with the BABAR detector at the PEP-II asymmetric- 
energy e + e~ collider at SLAC. We find evidence for B° — > p°p° with 3.5<r significance and measure 
the branching fraction B = (1.07 ± 0.33 ± 0.19) x 10 -6 and longitudinal polarization fraction ft = 
0.87 ± 0.13 ± 0.04, where the first uncertainty is statistical, and the second is systematic. The 
uncertainty on the CKM unitarity angle a due to penguin contributions in B — -> pp decays is 18° 
at the la level. We also set upper limits on the B° -» p°/ (980) and B° -» f (980) f (980) decay 
rates. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



Measurements of CP-violating asymmetries in the 
B°B° system test the flavor structure of the stan- 
dard model by over-constraining the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix [l[. The time- 
dependent CP asymmetry in the decays of B° or B° 
mesons to a CP eigenstate dominated by the tree- 
level amplitude b — > uud measures sin2a ff, where a c s 
differs from the CKM unitarity triangle angle a = 
arg [~V td V t * b /V ud V* b ] by a quantity Aa accounting for 
the contributions from loop (penguin) amplitudes. The 
value of Aa can be extracted from an analysis of the 
branching fractions of the B decays into the full set of 
isospin-related channels @. 

Branching fractions and time-dependent CP asymme- 
tries in B — > 7T7T, p7r, and pp have already provided in- 
formation on a. Since the tree contribution to the B° — > 
p°p° decay is color-suppressed, the decay rate is sensitive 
to the penguin amplitude. The B° 
much smaller branching fraction than B° 



p°p° decay has a 



P P 



and 



L £|, and therefore a 
7, Gil. This makes 



B+ -► p+p° channels [4J, i, 0, 0, 
stringent limit on Aa can be set 
the pp system particularly effective for measuring a. 

In B — > pp decays the final state is a superposition 
of CP-odd and CP-even states. An isospin-triangle re- 
lation 0] holds for each of the three helicity amplitudes, 
which can be separated through an angular analysis. The 
helicity angles 9\ and 9 2 are denned as the angles between 
the direction of ir + and the direction of the B in the rest 
system of each of the p° candidates. The resulting angu- 
lar distribution d 2 T /(T dcosOi dcos9 2 ) is 

| {jC 1 - h) sin 2 9i sin 2 9 2 + f L cos 2 6 X cos 2 9 2 \ , (1) 



2 ) is the longitudinal polariza- 
i,o. +i are the helicity amplitudes. 



where f L = \A \ 2 /(E\A X 
tion fraction and A\ 

In this paper we present the first evidence for the 
B° — > p°p° decay, the measurement of the longitudi- 
nal polarization fraction in this decay, and updated con- 
straints on the penguin contribution to the measurement 
of the unitarity angle a. 

These results are based on data collected with the 
BABAR detector |lll| at the PEP-II asymmetric-energy 
e + e~ collider [12]. A sample of 383.6 ± 4.2 million 



BB pairs was recorded at the T(4S) resonance with 
the center-of-mass (cm.) energy yfs = 10.58 GeV. 
Charged-particle momenta and trajectories are measured 
in a tracking system consisting of a five-layer double- 
sided silicon vertex tracker and a 40-layer drift chamber, 
both within a 1.5-T solenoidal magnetic field. Charged- 
particle identification is provided by measurements of the 
energy loss in the tracking devices and by a ring-imaging 
Cherenkov detector. 

We select B — > M\M 2 — > (7r + 7r~)(7r+7r~) candidates, 
with Mi t 2 standing for p° or /o candidate, from neu- 
tral combinations of four charged tracks that are con- 
sistent with originating from a single vertex near the 
e + e~ interaction point. We veto tracks that are pos- 
itively identified as kaons or electrons. The identifica- 
tion of signal B candidates is based on several kinematic 
variables. The beam-energy-substituted mass, tties = 
[(s/2 + pi-pB) 2 /Ef — p^] 1 / 2 , where the initial e + e~ four- 
momentum (Ei, pi) and the B momentum pe are denned 
in the laboratory frame, is centered near the B mass with 
a resolution of 2.6 MeV for signal candidates. The dif- 
ference AE — E™ — s/s/2 between the reconstructed 
B energy in the cm. frame and its known value sfsjl 
has a maximum near zero with a resolution of 20 MeV 
for signal events. Four other kinematic variables describe 
two possible 7r + 7r~ pairs: invariant masses mi, mi and 
helicity angles 6\, B 2 - 

The selection requirements for signal candidates are 
the following: 5.245 < m ES < 5.290 GeV, \AE\ < 
85 MeV, 550 < m%,2 < 1050 MeV, and |cos0i, 2 | < 0.98. 
The last requirement removes a region corresponding to 
low-momentum pions with low and more uncertain re- 
construction efficiency. In addition, we veto the copious 
decays B° -> D^-7r+ -> (h+ 7T 7r )-7T + , where h + refers 
to a pion or kaon, by requiring the invariant mass of the 
three-particle combination to differ from the D-meson 
mass by more than 13.2 MeV, or 40 MeV if one of the 
particles is consistent with a kaon hypothesis. 

We reject the dominant e + e~ — * qq (<j> = u, d, s, c) 
(continuum) background by requiring Icos^tI < 0.8, 
where 9t is the angle between the P-candidate thrust 
axis and that of the remaining tracks and neutral clusters 
in the event, calculated in the cm. frame. We further 
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suppress continuum background using a neural network 
discriminant £ , which combines a number of topological 
variables calculated in the cm. frame. Among those 
are the polar angles of the B momentum vector and the 
S-candidate thrust axis with respect to the beam axis. 
Other discriminating variables include the two Legcndre 
moments Lq and L2 of the energy flow around the B- 
candidate thrust axis 13j and the sum of the transverse 
momenta of all particles in the rest of the event, calcu- 
lated with respect to the B direction. 

After application of all selection criteria, iV canc i = 
64843 events are retained. On average, each se- 
lected event has 1.05 signal candidates, while in Monte 
Carlo 14] samples of longitudinally and transversely po- 
larized B° — > p°p° decays we find 1.15 and 1.03 candi- 
dates, respectively. When more than one candidate is 
present in the same event, the candidate having the best 
X 2 consistency with a single four-pion vertex is selected. 
Simulation shows that 18% of longitudinally and 4% of 
transversely polarized B° — > p°p° events are misrecon- 
structed with one or more tracks not originating from 
the B° — > p°p° decay. These are mostly due to combina- 
torial background from low-momentum tracks from the 
other B meson in the event. 

Further background separation is achieved by the use 
of multivariate £>-flavor-tagging algorithms trained to 
identify primary leptons, kaons, soft pions, and high- 
momentum charged particles from the other B [HI . The 
discrimination power arises from the difference between 
the tagging efficiencies for signal and background in seven 
tagging categories (c tag = 1..7). 

We use an unbinned extended maximum likelihood fit 



to extract the B° 



p°p° event yield and fraction of 



longitudinal polarization f^. We also fit for the event 
yields of B° — > p°f and B° — > f f decays, as well as of 
several background categories. The likelihood function is 



C exp n k II "■' ^J'(^i) 



(2) 



where rij is the unconstrained number of events for 



each event type j (B° -> p Q p° , B° -> p°/o(980), 
B° — ► /o (980) /o (980), three background components 
from B decays, and continuum), and Vj(xi) is the 
probability density function (PDF) of the variables 
xi = {toes, AE, £, ui\ , 777,2 1 cos 0i, cos 82, ct ag }i for the ith 
event. 

We use simulated events to parameterize the back- 
ground contributions from B decays. The charmless 
modes are grouped into two classes with similar kine- 
matic and topological properties: B° — > afir T and a 
combination of other charmless modes, including B° — > 
p°K*°, B+ -> p+p Q , B -► pir, and B° -► p+p~. One 
additional class accounts for the remaining neutral and 
charged B decays to charm modes. We ignore any other 



four-pion final states whose contributions are expected 
to be small in our invariant mass window. 

Since the statistical correlations among the variables 
are found to be small, we take each Vj as the product 
of the PDFs for the separate variables. Exceptions are 
the kinematic correlation between the two helicity angles 
in signal, and mass-helicity correlations in other i3-decay 
classes and misreconstructed signal. 

We use double-Gaussian functions to parameterize the 
toes and AE PDFs for signal, and a relativistic Breit- 
Wigner functions for the resonance masses of p° and 
/o(980) 0. The angular distribution at production 
for B° -> p°p°, B° -► p°f , and B° -► / / modes 
(expressed as a function of the longitudinal polariza- 
tion in Eq. ([1]) for B° — > p°p°) is multiplied by a de- 
tector acceptance function C/(cos6>i, cos6* 2 ), determined 
from Monte Carlo. The distributions of misreconstructed 
signal events are parameterized with empirical shapes in 
a way similar to that used for B background discussed 
below. The neural network discriminant £ is described 
by three asymmetric Gaussian functions with different 
parameters for signal and background distributions. 

The PDFs for non-signal B decay modes are generally 
modeled with empirical analytical distributions. Several 
variables have distributions identical to those for signal, 
such as ttt-es when all four tracks come from the same B, 
or 7r + 7r _ invariant mass mi,2 when both tracks come from 
a p° meson. Also for some of the modes the two ir + ir~ 
pairs can have different mass and helicity distributions, 
e.g. when only one of the two combinations comes from 
a genuine p° or /o meson, or when one of the two pairs 
contains a high-momentum pion (as in B — ► ai7r). In 
such cases, we use a four- variable correlated mass-helicity 
PDF. 

The signal and B-background PDF parameters are ex- 
tracted from simulation. The Monte Carlo parameters 
for toes, AE, and £ PDFs are adjusted by compar- 
ing data and simulation in control channels with simi- 
lar kinematics and topology, such as B° — > D~ir + with 
D~ — > K + n~Tt~. The continuum background PDF pa- 
rameters are left free in the fit. Finally, the discrete B- 
fiavor tagging PDFs for signal modes are obtained in ded- 
icated fits to events with identified exclusive B decays. 
The tagging PDFs for inclusive B backgrounds are deter- 
mined by Monte Carlo and their systematic uncertainties 
are studied in data. 

Table U shows the results of the fit. The B° -> p a p° 
decay is observed with a significance of 3.5er, as deter- 
mined by the quantity y— 2 log(£o/£max), where £ max 
is the maximum likelihood value, and £0 is the likelihood 
for a fit with the signal contribution set to zero. It cor- 
responds to a probability of background fluctuation to 
the observed signal yield of 2 x 10 -4 , including system- 
atic uncertainties, which are assumed to be Gaussian- 
distributed. We do not observe significant event yields 
for B° -> p /o(980) and B° -> / (980)/ (980) decays. 
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FIG. 1: Projections of the multidimensional fit onto (a) jties, 
(b) AE, (c) di-pion invariant mass (mi is shown, distribution 
of 77i2 is similar), and (d) cosine of the helicity angle (cos Si 
is shown), after a requirement on the signal-to-background 
probability ratio with the plotted variable excluded. This re- 
quirement enhances the fraction of signal events in the sample. 
The data points are overlaid by the solid black line, which cor- 
responds the full PDF projection. The individual B° — -> p a p° 
PDF component is also shown with a solid red line. The sum 
of all other PDFs (including B° -> p°f and B° -> f fo com- 
ponents) is shown as the dashed blue line. The D-meson veto 
causes the acceptance dip seen in (d). 



Background yields are found to be consistent with ex- 
pectations. In Fig. [1] we show the projections of the fit 
results onto ttieSi AiS, mi, and cos#i variables. 

Dominant systematic uncertainties in the fit originate 
from statistical errors in the PDF parameterizations, due 
to the limited number of events in the control samples. 
The PDF parameters are varied by their respective un- 
certainties to derive the corresponding systematic errors 
(±10, lg, ±4 events for p°p°, p°fo, and f f respectively, 
and 0.03 for fa). We also assign a systematic error of 2 
events for p°p°, 3 events for p°fo, and 1 event for /o/o 
(0.01 for fi) to account for a possible fit bias, evaluated 
with Monte Carlo experiments. The above systematic 
uncertainties do not scale with event yield and are in- 
cluded in the calculation of the significance of the result. 

We estimate the systematic uncertainty due to the in- 
terference between the B° — ► p Q p a and af ir T decays us- 
ing simulated samples in which the decay amplitudes for 
3° — > p°p are § 
and those for B° — 
fraction of (33.2 ±4. 



B u — > p°p are generated according to this measurement 



correspond to a branching 
6 ITJ. Their amplitudes are 
modeled with a Breit-Wigner function for all p — > tttt 



x 10" 



TABLE I: Summary of results: event yields (n); fraction of 
longitudinal polarization (/l); selection efficiency (Eff) corre- 
sponding to measured polarization; branching fraction (23 s ig); 
branching fraction upper limit (UL) at 90% CL; and signifi- 
cance including systematic uncertainties. The systematic er- 
rors are quoted last. We also show the background event 
yields for ai7r, qq, charmless, and other BB components (sta- 
tistical uncertainties only). 



Quantity 



Value 



n(B° - p°p°) 

h 

Eff (%) 
& ig (xlO" 6 ) 

Significance, stat. only (<r) 
Significance, syst. included (a) 



100 ± 32 ± 17 
0.87 ±0.13 ±0.04 

24.2 ± 1.0 
1.07 ±0.33 ±0.19 
3.7 
3.5 



n(B° -» p°f ) 
Eff (%) 

B Bi g XB(f ^TT+TT") (XlO" 6 ) 
ULxB(/ -> 7T+7T-) (xlO" 6 ) 



+7 



20 ±21 
26.1 ± 1.6" 
0.19 ±0.21 ta°l 
0.53 



77(i3 -» / /o) 

Eff (%) 

S sig xB 2 (f ^TT+TT-) (XlO" 6 ) 
ULx/3 2 (/o -» TT + tO (XlO" 6 ) 



-3±9±5 
28.6 ± 1.1 
-0.03 ±0.08 ±0.04 
0.16 



n{B° a+TT*) 
n(charmless) 
n(BB) 
n(qq) 



81 ±25 

-i 7 +107 
11 -96 

3198 ± 224 
61469 ±311 



and a\ —> pit combinations and their relative phase is as- 
sumed to be constant across the phase space. The strong 
phases and CP content of the interfering state a^Tr^ are 
varied between zero and a maximum value using uniform 
prior distributions. We take the RMS variation of the 
average signal yield (14 events for the p°p° yield, or 0.03 
for /i) as a systematic uncertainty. 

Uncertainties in the reconstruction efficiency arise from 
track finding (2%), particle identification (2%), and other 
selection requirements, such as vertex probability (2%), 
track multiplicity (1%), and thrust angle (1%). 

To constrain the penguin contributions to B — > pp 
decays, we perform an isospin analysis, by minimizing 
a x 2 term that includes the measured quantities ex- 
pressed as the lengths of the sides of the isospin trian- 
gles. We use the measured branching fractions and frac- 
tions of longitudinal polarization of the B + 



and B 



o 



P + P° 



P P~ 



[jj decays, the CP-violating parameters 
S~^~ and ~ determined from the time evolution of the 
longitudinally polarized B° — > p + p~ decay Q, and the 
branching fraction and polarization of B° — > p°p° from 
this analysis. We assume uncertainties to be Gaussian 
and neglect 1=1 isospin contributions, electroweak loop 
amplitudes, non-resonant and isospin-breaking effects. 

With the B° — > p°p° measurement we obtain a 68% 
(90%) CL limit on |Aa| = |a - a off | < 18° (< 20°). 



7 




FIG. 2: Ax 2 as a function of Aa obtained from the isospin 
analysis discussed in the text. The dashed lines at A\ 2 = 1 
and Ax 2 = 2.7 are taken for the la (68%) and 1.64a (90%) 
interval estimates. 



Fig. [2] shows A\ 2 as a function of Aa. The cen- 
tral value of a obtained from the isospin analysis is 
the same as a c g, which is constrained by the relation 
sin(2a c ff) = S^~/(l - C^ -2 ) 1 / 2 and is measured with 
the B° -> p+p- decay 0. 

The error due to the penguin contribution becomes 
the dominant uncertainty in the measurement of a using 
B — > pp decays. However, once the sample of B° — > p°p° 
decays becomes more significant, time-dependent angu- 
lar analysis will allow us to measure the CP parameters 
5°° and C°°, analogous to and C^~, resolving am- 
biguities inherent to isospin triangle orientations. 

In summary, we find evidence for B° — > p°p° decay 
with 3.5ct significance. We measure the B° p°p° 
branching fraction of (1.07±0.33±0.19) x 1(T 6 and deter- 
mine the longitudinal polarization fraction for these de- 
cays of /l = 0.87±0.13±0.04. The measurement of this 
branching fraction combined with that for B + — * P°P + 
and B° — ► p + p~ decays provides a constraint on the 
penguin uncertainty in the determination of the CKM 
unitarity angle a. These results supersede our previous 
measurements 0] ■ We find no significant evidence for the 
decays B° -► p°f and B° -> / / . 
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